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Abstract. The dust extinction towards bright H ii region- 
sin NGC 598 and NGC 5457 has been studied in detail by 
forming line ratios of Balmer and Paschen emission lines 
covering a large wavelength range. Three homogeneous 
models of the geometrical distribution of the emitting 
sources and the obscuring dust have been tested. Only for 
low extinctions can the data be fit by a homogeneous slab. 
For most of the observed H ii regions the Witt et al. (1992) 
’dusty nucleus’ model matches the observations equally 
well as the usual assumption of a foreground screen, but 
the former implies much larger actual dust contents. Spa¬ 
tial variations in the optical depth in V of the order 0.3 
across a region are found. 

Key words: galaxies: spiral - galaxies: ISM - dust, ex¬ 
tinction - HII regions 


1. Introduction 

Investigations of external galaxies are often faced with the 
problem of disentangling the effects of dust and the phys¬ 
ical properties, e.g. age, IMF, chemical composition etc. 
of the emitting source. 

The effect of dust extinction goes roughly as and 
will consequently change the observed spectral energy dis¬ 
tribution. Studies and interpretation of the physical prop¬ 
erties in external galaxies therefore rely on how well we 
can account for the attenuation of the emitted light due 
to dust. The wavelength dependence makes multiwave¬ 
length analyses favourable and a growing number of in¬ 
vestigations combining optical and IR wavelengths have 
been made both using imaging (e.g. Jansen et al. 1994; 
Witt et al. 1994; Evans 1995; van Driel et al. 1995) and 
spectroscopy (e.g. Puxley & Brand 1994; Galzetti et al. 
1996). 

Send offprint requests to\ larsp@obs.aau.dk 

* Based on observations collected at the Nordic Optical Tele¬ 
scope, La Palma, Spain 


In the analysis of spectra of extragalactic Hii re¬ 
gions the correction for dust extinction is usually made on 
the assumption that the intrinsic hydrogen emission line 
spectrum is known through case B recombination and any 
difference between the emergent and intrinsic flux ratios 
can be ascribed to dust obscuration following the Galactic 
extinction curve. 

Throughout the years, the use of the ratio of the two 
strongest Balmer lines, Ha/H/3, sometimes in combination 
with Hy, has been prevailing, albeit they only span a lim¬ 
ited wavelength range, because they will be included in 
almost any optical spectra and are easy to detect. It has 
long been proposed to expand the baseline for extinction 
determination by comparing near-IR hydrogen Paschen 
lines with emission lines from the Balmer series at short 
wavelengths (cf. Greve et al. 1989; Osterbrock 1989). Such 
line pairs are all separated by wide wavelength intervals 
over which the extinction has a large effect and is therefore 
potentially easier to deduce. Additionally, some line pairs 
form corresponding multiplet lines, Pn/Hn originating at 
the same upper atomic level with relative strength depend¬ 
ing primarily on the transition probability and minimised 
dependence on theoretical recombination line calculations 
(Greve et al. 1994). 

In Petersen & Gammelgaard (1996; hereafter Paper 
I) we have demonstrated the feasibility of observing lines 
from the Paschen series and the short end of the Balmer 
series from giant extragalactic H ii regions simultaneously 
with a two-channel spectrograph. In the present paper we 
will further discuss dust extinction towards H ii regions in 
NGG 598 and NGG 5457 as derived from spectra obtained 
with the same instrument. 

The straightforward application of the Galactic inter¬ 
stellar extinction curve derived from point sources to ex¬ 
tended emission regions in external galaxies is, however, 
problematic, because this procedure implicitly assumes a 
dust configuration of a homogeneous foreground absorb¬ 
ing screen. Intuitively one would expect things to be more 
complicated in reality with dust distributed in a patchy or 
clumpy way maybe mixed with the emitting source (Witt 
et al. 1992; hereafter WTG; Galzetti et al. 1994, 1996; 
Puxley & Brand 1994) and with scattering of, particu- 
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larly blue, photons into the line of sight by dust grains 
(Bruzual et al. 1988; WTC). This could all lead to a much 
greyer effective extinction curve than the standard Galac¬ 
tic curve. This was exemplified in Paper I by using an 
extinction curve with Ry = Ay/Eb-v = 6, here we will 
address this issue by fitting three different models of rel¬ 
ative spatial distribution of emitting source and dust to 
the measured line ratios. 


produce less reddening in a given geometry with an ex¬ 
tended source distribution and approach lower finite lim¬ 
iting reddening with increasing dust mass (Witt & Gordon 
1996). Thus the dust content derived on the assumption 
of homogeneous distributions can underestimate the true 
dust mass. 

2.1. Foreground screen 


2. The models 

When data from a wide wavelength range are available it 
is possible to test various models for the geometrical dis¬ 
tribution of the dust. Although the data presented here 
only reach PJ at 1.0 fxm and in two cases Py at 1.1 /rm, 
they can still be used to test some basic models and have 
the advantage of being recorded simultaneously with the 
same instrument as the Balmer emission lines. Three mod¬ 
els have been considered: (1) The naive but widely used 
uniform foreground extinction screen, (2) a slab model of 
homogeneous mixture of dust and emitting sources and 
(3) that of the WTC scenarios which best models an ex- 
tragalactic Hii region, the so-called ’’dusty nucleus”. 

The extinction can be quantified in many ways, e.g. 
as the logarithmic extinction at H/3, C(RP), the visual 
extinction in magnitudes, Ay or as the optical depth in 
V, Ty, depending on the purpose it is used for. In this 
work we will express the derived dust extinction as Ty for 
uniformity among the models. In the case of a foreground 
screen simple relations exist between these, with Ty = 
1.98C'(H/3) = Av/1.086. 

Depending on the actual distribution the same total 
amount of dust can produce very different attenuation of 
the emitted light. Therefore we define an effective optical 
depth Teff = — hi{lFohs/IFem)■ In most cases the actual op¬ 
tical depth will exceed Tes, thus some source distributions 
are capable of hiding away larges amount of extinguishing 
matter. 

More complicated models with dumpiness could be in¬ 
cluded, but due to the lack of near-IR lines at longer wave¬ 
lengths and the rather low extinction for some of the re¬ 
gions we do not find it justified to attempt to discriminate 
between clumpy and homogeneous distributions. Thus we 
limit ourselves to the three homogeneous models but it 
must be stressed that these will produce lower limits to the 
dust content as compared to corresponding clumpy mod¬ 
els. Calculations of the effect on the emitted radiation field 
by a clumpy screen with dust-free interclump medium ex¬ 
cluding scattering (Natta & Panagia 1984), a two-phase 
slab (Boisse 1990) and spherical two-phase clumpy sys¬ 
tems with non-conservative scattering (Witt & Gordon 
1996) all show reduced overall effective optical depth rel¬ 
ative to homogeneous distributions of the same mass and 
flatter effective extinctions curves. As clumps will become 
optically thick much quicker than an equivalent homoge¬ 
neous distribution would, a clumpy dust distribution will 


This is the most simple geometry of dust distribution, in¬ 
herited from stellar astronomy but likely to be invalid for 
many extragalactic studies where more complicated ge¬ 
ometries are expected (WTC). On the other hand Calzetti 
et al. (1996) found evidence that it might be applicable for 
star burst environments, where large fractions of the dust 
have been depleted or ejected. By chance it happens to be 
the geometry where any amount of dust has the strongest 
effect on the emitted flux and taken at face value it will 
always yield a lower limit to the true dust content. 

For any pair of emission lines the relationship of ob¬ 
served line flux ratio and the predicted ratio in the absence 
of dust is 

^^grx2-rxy (1) 

-Tip 

Hence, as an exception, the foreground screen has by defi¬ 
nition T = Teff. Introducing a normalised extinction curve 
TxjTy = A\/Ay one can solve for Ty 


Ty 


h\{Ro/Rp) 
TXljTy — TAl/rv 


2.51og(A’o/Rp) 

AxijAy — AwjAy ' 


( 2 ) 


if transforming to magnitude. The first expression in Eq. 
(2) is of cause the motivation for illustrating the data as 
in Fig. 2-4. 


2.2. Homogeneous slab 

A slab model of mixture of dust and emitting sources may 
resemble the actual geometry better than the overlying 
screen in some cases, although a single or a few knots of 
emission surrounded by non-uniform dust is probably the 
dominating picture. For a homogeneous mixture through¬ 
out neglecting scattering Ro/Rp can be described in terms 
of the optical depth (Disney et al. 1982, Puxley & Brand 
1994) 


^ ^ / 1 - e \ 

Rp tai V1 - ) 


( 3 ) 


When introducing relative values TxjTy we end up with 
the only free parameter Ty, which is fit to give the best 
match to the data points. 
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2.3. Dusty nucleus 

This is probably the most realistic model of the real dust 
distribution affecting the emitted flux from an H ii region. 
The WTC ’dusty galactic nucleus’ environment consists 
of a dust-free spherical source surrounded by a dust-shell 
without emitting sources and is characterised by the op¬ 
tical thickness of the shell, ry. The main difference be¬ 
tween this model and the other two is the presence of a 
dust-shell subtending a large solid angle as seen by the 
source and the inclusion of scattered light back into the 
line of sight from it, whereas the other models only in¬ 
clude pure extinction (absorption and scattering out of 
the line of sight). This difference has increasing impact 
in the short-wavelength end resulting in a ’bluing’ of the 
spectrum partly cancelling the reddening effect of the ex¬ 
tinction and leading to a flatter effective extinction curve. 
Consequently, the observational effect will be to see values 
of ln(i?o/.Rp) close to zero in the blue end of Figs. 2-4 even 
at large optical depths in contrast to the two other distri¬ 
butions, which lead to decreasing values of ln(i?o/i?p) as 
function of Txfry with a slope depending on ry. 

WTC tabulate the amount of direct and scattered 
light received at the central wavelength of the broad band 
colours calculated from Monto Carlo simulations of the 
transfer of radiation. In this analysis we have interpolated 
their tables to the wavelength of the relevant emission lines 
by cubic splines and use the sum of the two contributions 
as the observed flux. 


3. Observations 

The spectra were obtained with the Arhus-Tromsp Low 
Dispersion Spectrograph at the Nordic Optical Telescope, 
La Palma on July 31 and August 1 1994, together with 
some of the data presented in Paper 1. 

Three extended H ii regions in NGC 598 all classified 
as extremely strong by Boulesteix et al. (1974) were ob¬ 
served. NGC 604 is the largest Hii region in NGC 598 
and resembles in some respect the 30 Dor complex in the 
LMC. The elemental abundance has been studied in de¬ 
tail by Diaz et al. (1987) using multiple slit positions. We 
have positioned the slit in an E-W direction through the 
center of the complex and extracted spectra from three 
bright knots roughly coinciding with their areas denoted 
D, B and C. In NGC 595, another large region, the slit 
was positioned across the southern part. B88 shows a fee¬ 
ble ring-like structure and one spectrum on the east and 
west side of the ring was extracted. 

The five Hii regions observed in NGC 5457 are all sit¬ 
uated in the south-east part of the galactic disk. The giant 
complex NGC 5461 was observed in two slit positions both 
centered visually on the peak luminosity from short Ha- 
exposures. One position (Jul 31) is E-W while the other 
is roughly in the NE-SW direction to include H1026 and 


Table 1. Journal of observations 


Date 

Galaxy 

Hii reg.’" 

Apertures 

(") 

Total exp. 

(s) 

Aug 1 

NGC 598 

N595 

2.2 

X 33.4 

3600 

Aug 1 

NGC 598 

B88 

2.2 

X 8.3 

3600 




2.2 

X 8.3 

3600 

Aug 1 

NGC 598 

N604 

2.2 

X 9.2 

1800 




2.2 

X 7.5 

1800 




2.2 

X 8.3 

1800 

Jul 31 

NGC 5457 

N5461 

2.2 

X 15.8 

3600 

Jul 31 

NGC 5457 

H949 

2.2 

X 10.8 

1800 



H961 

2.2 

X 8.3 

1800 

Aug 1 

NGC 5457 

H1026 

2.2 

X 8.3 

6000 



N5461 

2.2 

X 10.0 

6000 



H1117 

2.2 

X 9.2 

6000 


^ ID numbers from Boulesteix et al. (1974) and Hodge et 
al. (1990) are preceded by B and H, respectively 


H1117 with NGC5461. A complete log of the observations 
is given in Table 1. 

During both nights the spectrophotometric standard 
star BD -|-17°4708 (Oke & Gunn 1983) was observed mul¬ 
tiple times. For all Hii regions and the standard star the 
blue and near-IR spectral ranges were observed simulta¬ 
neously in the blue and red channel of the spectrograph to 
ensure precise relative spectrophotometry of the Balmer 
and Paschen lines. The two channels have dispersions of 
2.0 A/pixel and 7.5 A/pixel and cover 3950-4900 A and 
8000-11000 A, respectively. A detailed description of the 
instrument and observation and reduction procedures can 
be found in Paper 1. 

3.1. Data reduction 

All the CCD exposures were reduced with standard IRAF 
packages for flat fielding, calibrations and extraction of 
the one-dimensional spectra. Rectification of the two- 
dimensional frames and transformation to a common spa¬ 
tial coordinate system were done with a special slit mask 
as explained in Paper 1. 

In the combined object frames extraction apertures 
were set to enclose each distinct emission region; the re¬ 
sulting apertures are listed in Table 1. The sky background 
was sampled along large part of the slit, avoiding any sign 
of emission, to get a good S/N of the sky subtraction in 
the near-IR. 

3.2. Relative line fluxes 

In Fig. 1 we present an example of the final flux calibrated 
spectra in the blue and near-IR wavelength range from of 
one of the H ii regions with high signal. The line fluxes 
have been measured by fitting Gaussian profiles to the line 
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Fig. 1. Left panel, Blue spectrum of NGC 5461 showing four strong Balmer emission lines. Right panel, near-IR spectrum of 
the same region. The Paschen lines PI5-P7 are clearly identihed, although with rapidly increasing noise for A > 10500 A. Note 
the different scales on the axes 


profiles or two Gaussians when deblending He from the 
nearby [Ne ill] line. As evident from Fig. 1 the spectrum 
has low S/N at the position of P 7 . This line was only 
detected in two of the spectra and the flux could not in 
any case be determined from the IRAF procedure. Instead 
the flux was measured by fitting a Gaussian with fixed Ag. 

For the strong Balmer lines we estimate the accuracy 
of the flux measurements to be of the order 5 %. Because 
He is measured from a deblending procedure it has a larger 
flux of error estimated to 20 %. In the near-IR flux uncer¬ 
tainties are dominated by the more critical sky subtrac¬ 
tion and determination of the spectral continuum level. 
At P15 we tend to reach the series limit where the faint 
Paschen lines start to blend giving rise to a false con¬ 
tinuum, this line is therefore excluded from the present 
analysis. Judged upon the line profile fitting procedure we 
assign an error estimate of 20 % to the two faintest lines 
included (P14 and P13) and to P9, which is strongly af¬ 
fected by the atmospheric H 2 O absorption band, and half 
of this to P12-P10 and PJ. 

The nebular Balmer emission lines are affected by the 
absorption lines in the underlying stellar spectrum. This 
will suppress the observed flux of the lines and relative 
more so at the short end of the line series, thus partly 
counteracting the effect of scattering on dust grains. The 
exact correction for underlying Balmer absorption is very 


uncertain, because it depends on the actual stellar popu¬ 
lation and age. Due to the lack of any detailed knowledge 
we have decided to use the average absorption equiva¬ 
lent width of 1.9 A, found by McGall et al. (1985), for 
all Balmer lines. A similar value was found by Diaz et al. 
(1987) in their analysis of NGC 604. In most cases the 
correction is quite small and will be less than the quoted 
error bars, except for HJ in the faintest regions and at He, 
where it may constitute a fair fraction of the measured 
line flux. We deal with this problem by assigning errors to 
include no correction and corrections up to 4 A equivalent 
width in a fashion similar to Puxley & Brand (1994). The 
large uncertainty for some of the spectra further cautions 
the usage of He. 

In Table 2 we list the line flux relative to H/3 of all 
identified hydrogen emission lines of the observed H ii re¬ 
gions in NGC 598 and NGC 5457. 


4. Discussion 

The quotient of observed and predicted line fluxes in three 
of the observed regions relative to H/3, 


Ro [ -^(A) \ /[ /(A) \ 

Rp {HmJo/ vimJp 


( 4 ) 
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Table 2. Observed line fluxes (normalised at = 1) 


ID 






Ao = ( 

T(A) \ 

/ o 







N595 

B88,l 

B88,2 

N604,l 

N604,2 

N604,3 

N5461 

H949 

H961 

H1026 

N5461 

H1117 

He 

0.102 

0.130 

0.122 

0.137 

0.125 

0.096 

0.104 

0.122 

0.098 

0.091 

0.101 

0.188 

H5 

0.232 

0.235 


0.253 

0.256 

0.252 

0.209 

0.210 

0.212 

0.178 

0.191 


Hy 

0.409 

0.437 

0.491 

0.441 

0.456 

0.451 

0.392 

0.402 

0.395 

0.387 

0.365 

0.370 

Hd 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

P14 

0.011 



0.010 

0.009 


0.016 






P13 

0.019 

0.009 


0.013 

0.012 


0.020 




0.022 


P12 

0.023 

0.014 

0.013 

0.014 

0.015 


0.020 




0.034 


Pll 

0.029 

0.015 

0.020 

0.017 

0.018 


0.025 




0.039 


PIO 

0.048 

0.022 

0.044 

0.021 

0.021 

0.029 

0.035 


0.049 

0.062 

0.039 

0.115 

P9 

0.064 

0.034 

0.052 

0.038 

0.037 

0.055 

0.059 

0.057 

0.068 

0.110 

0.071 

0.126 

P5 

0.159 

0.085 

0.164 

0.090 

0.101 

0.064 

0.158 

0.145 

0.200 

0.360 

0.211 


Py 





0.378 


0.297 






J^(H/3)“ 

-12.55 

-13.07 

-13.48 

-12.43 

-12.61 

-13.59 

-12.73 

-13.38 

-13.50 

-13.90 

-12.79 

-14.43 


^ Logarithm of the H/f-flux in units of erg cm ^ s ^ 


are plotted in Figs. 2-4 with ta/tv = Ax/Ay on the ab¬ 
scissa as derived from the fitting formula by Cardelli et al. 
(1989) with the standard value Ry = 3.1. The choice of 
normalising to H/3 is simply a matter of convenience and 
when deriving ry from the corresponding multiplet lines 
the relevant Balmer lines are used in the general expres¬ 
sion Eq. (2). 

The theoretical line ratios Rp for case B recombination 
are adopted from Hummer & Storey (1987). The electron 
temperature and density have been determined by Kwit- 
ter & Aller (1981), Rayo et al. (1982), Diaz et al. (1987) 
and Vilchez et al. (1988). For the two regions (H1026 and 
H1117) where Tg and Ug can not be acquired from the lit¬ 
erature we used typical values of 10"^ K and 10^ cm“^ and 
include realistic variations in the error estimates in Tables 
3 and 4. The weak variation of Rp with these parameters 
will in any case only introduce uncertainties smaller than 
the expected observational errors (Greve et al. 1994). 

According to Eq. (2) the overlying screen will manifest 
itself as a straight line with a slope equal to —ry in this 
type of plot. Figure 2 illustrates this for NGC 5461. The 
best linear fit to the data points gives ry = 1.04 (Table 
3) with = 13, but fails to match the data points at the 
longest wavelengths, i.e. smallest Txjry. More complicated 
geometries introduce curvature in the plots (see Eq. (3) 
and Figs. 2-3). A homogeneous slab with optical depth 
Ty = 4.3 fits the data equally well or perhaps even better 
(x^ = 11), especially at the long wavelength end. The 
’dusty nucleus’ model of ry = 3 does not fit quite as well 
(x^ = 34). Although it gives a good match to Py and 
P(5 it exaggerates somewhat the ’blueing’ of the spectrum 
at short wavelengths caused by scattering. Bearing the 


problems of the He line flux in mind this might not be 
significant, though. 

Given the indicated uncertainties it would be wrong 
to rule out any of the analysed dust configurations for 
this particular H ii region, but it is clearly demonstrated 
that while the more complex dust geometries have effec¬ 
tive optical depths only moderately higher than a fore¬ 
ground screen the true dust content can be as much 3-4 
times higher. A definite discrimination of various models 
will need observations of emission lines in the Paschen and 
Brackett series at still longer near-IR wavelengths (Puxley 
& Brand 1994). 

From the observation of NGC 5461 in the other slit 
position the following night we derive larger values of the 
optical depth, ry = 1.41 for the foreground screen and 
Ty = 5 for the ’dusty nucleus’. No acceptable fit was pro¬ 
duced by the homogeneous slab model. This deviation is 
larger than the estimated uncertainties and we believe it 
is real and a consequence of the non-coinciding slit posi¬ 
tions, and only partly attributable to the effect of differ¬ 
ential refraction (Filippenko 1982) and the overall lower 
signal, since the principal Paschen lines all stand well out 
of the noise level of the continuum and the slit is much 
larger than the typical seeing encountered during the ob¬ 
servations (< 1") and was close to the optimal orientation. 
As discussed in section 4.1 variations of this order of mag¬ 
nitude are likely to be present across giant H ii regions at 
the spatial scale sampled here. 

The optical depths for all the H ii regions observed in 
NGC 5457 deduced from the three different dust distri¬ 
butions are collected in Table 3. H949 is the only other 
region where it was possible to obtain a best fitting model 
from the homogeneous slab distribution with Ty < 10. 
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Fig. 2. Observed and predicted 
line flux relative to H/3 in the H ll 
region NGC5461 with three best 
fit models for the dust distribu¬ 
tion: Foreground screen (dotted 
line), homogeneous slab (dashed 
line) and the WTG ’dusty nucleus’ 
(dot-dashed) 



rx/rv - Aj/Ay 


Fig. 3. Same as Fig. 2 but for the 
region H1026 where the quotients 
of observed and predicted line ra¬ 
tios clearly cannot be fit by the 
homogeneous slab model 


(For larger values of ry the slab model have an asymp¬ 
totic behaviour which will not change the shape of the 
curve with increasing optical depth.) The data from this 
HII region are also well fit with a ’dusty nucleus’ model 
of ry = 3. For H1026 and H1117 we find remarkable high 
values of the optical depth that obviously cannot be fit by 
the homogeneous slab (see Fig. 3), and is best described 
by a WTC model with ry equal 4 or 5, suggesting that 
quite large amounts of dust are associated with the spiral 
arm along which these two regions are positioned. 

In Fig. 4 we plot the relative line fluxes of NGC 604 
in NGC 598 in the same fashion as for NGC 5461. Since 
Py is just barely identified it has not been included in 
the diagram. Evidently, this H ii region experiences much 
less dust extinction making it even harder to distinguish 
between any of the three models. The plotted curves fit 
the data points equally well with both the homogeneous 


slab model and the ’dusty nucleus’ model having ry = 
1, but apparently the scattering properties of the WTC 
model follow better the flattening of \n{Ro/Rp) for the 
blue optical lines. None of the models adequately fit the 
measurement of P(5, which indicate higher extinction. 

Table 4 lists the values of the optical depth derived 
from the three models for the observed H ii regions in 
NGC 598. In general they show low values of ry with the 
highest found for NGC 595. For this Hii region the obser¬ 
vations are best described by a homogeneous slab or the 
’dusty nucleus’ model again with indications of a ’blue¬ 
ing’ of the spectrum at short wavelengths. For the other 
observations none of the considered models stand out. 

Two pairs of corresponding multiplet lines, Py/HJ and 
P5/He, are included in the present observations but the 
former only for two of the Hii regions. While the usage 
of any of the two pairs are not unproblematic and involve 





L. Petersen et al.: Dust extinction of Hll regions in NGC 598 and NGC 5457 


7 


Table 3. Optical depth rv for the observed Hll regions in NGC 5457 , a indicate that no acceptable fit conld be obtained. 
The Galactic extinction is 0. in this direction (de Vauconleurs et al. 1991) 



N5461 

H949 

H961 

H1026 

N5461 

H1117 

Foreground screen 

Py/Hd 

1.27 ±0.2 






Pd/He 

1.35 ±0.2 

1.14 ±0.3 

1.65 ±0.2 

2.27 ±0.3 

1.65 ±0.2 


Pd/H/3 

1.35 ±0.2 

1.26 ±0.2 

1.68 ±0.2 

2.45 ±0.2 

1.73 ±0.2 


linear fit 

1.04 ±0.07 

1.15 ±0.06 

1.31 ±0.10 

1.96 ±0.2 

1.44 ±0.09 

2.26 ±0.3 

Homogenous slab 

4.3 

5.5 

- 

- 

- 

- 

Homogenous slab, eff 

1.5 

1.7 

- 

- 

- 

- 

WTC 

3 

3 

3 

4 

3 

5 

WTC, eff 

1.7 

1.7 

1.7 

2.4 

1.7 

3.1 

Peimbert & Spinrad (1970) 

0.40 ±0.1 




0.40 ±0.10 


Smith (1975) 

0.83 




0.83 


Rayo et al. (1982) 

1.49 

( 1.39 ) 


1.49 


McGall et al. (1985) 

1.16 ±0.13 




1.16 ±0.13 


Kennicutt & Garnett (1996) 

0.71 ±0.10 

0.77 ±0.06 


0.22 ± 0.04 

0.71 ±0.10 


Skillman & Israel (1988)^ 

1.23 




1.23 


Rosa & Benvenuti (1994)*^ 

0.87 




0.87 



Derived from the Br7/H/1 ratio 

Based on HST observations and population models 



rx/rv - Ax/Ay 


Fig. 4. As Fig. 2 but for the cen¬ 
tral emission region of the NGC 
604 complex 


large flux uncertainties they do on the other hand provide 
the advantage of common upper levels and less reliance on 
the recombination models. Here they have been included 
in Tables 3 and 4 for reason of comparison with deriva¬ 
tions based on all observed emission lines. It is evident 
that ratios of these Paschen and Balmer lines always yield 
a higher value of ry, although with larger errors than us¬ 
ing the full set of observed lines. This could be due to an 
underestimation of the He flux in the deblending process 
and the correction for underlying Balmer absorption or 
the low S/N at the position of Py reflected in the quoted 


error bars. But Figs. 2 and 4 also show that Py and Pd al¬ 
ways lie above the straight line outlined by the remaining 
points, thereby indicating some curvature away from the 
foreground screen model. This higher optical depth deter¬ 
mined from emission lines with long wavelength baselines 
is further discussed in section 4.2. 

Finally, we have included extinction calculations based 
on the single line pair made up of the Paschen and Balmer 
line with highest S/N, namely Pd/H/3. The derived opti¬ 
cal depths from this ratio agrees well with those found 
on the basis of the corresponding line pair Pd/He, expect 
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Table 4. Optical depth ry for the observed Hll regions in NGC 598. rv“* = 0.13 (de Vauconlenrs et al. 1991) 



N595 

B88,l B88,2 

N604,l 

N604,2 

N604,3 

Foreground screen 






Py/Hii 




1.31 ±0.4 


P5/He 

1.38 ±0.2 

0.54 ±0.4 1.22 ±0.3 

0.57 ±0.2 

0.76 ±0.2 

0.59 ±0.3 

P5/H/3 

1.35 ±0.2 

0.52 ±0.2 1.38 ±0.2 

0.62 ±0.2 

0.76 ±0.2 

0.18 ±0.2 

lin. fit 

1.08 ±0.08 

0.32 ±0.06 0.87 ±0.15 

0.35 ±0.06 

0.39 ±0.09 

0.41 ±0.14 

Homogenous slab 

6.5 

0.8 3.8 

0.9 

1.0 

1.0 

Homogenous slab, eff 

1.9 

0.4 1.4 

0.4 

0.5 

0.5 

WTC 

3 

1 2 

1 

1 

1 

WTC, eff 

1.7 

0.5 1.1 

0.5 

0.5 

0.5 

Peimbert & Spinrad (1970) 



( 

0.44 ±0.10 

) 

Smith (1975) 

1.00 


( 

0. 

) 

French (1980) 



( 

0.80 

) 

Kwitter & Aller (1981) 

1.35 

( 0.50 ) 

( 

0.75 

) 

McCall et al. (1985) 

0.95 





Diaz et al. (1987) 



0.59 

0.59 

0.79 

Melnick et al. (1987) 

0.77 


( 

0.59 

) 

Vilchez et al. (1988) 


( 0.20 ) 

( 

0.69 

) 


in the case of N604,3 due to the very small He flux mea¬ 
sured for this region. For the faintest regions P5/H/3 gives 
a smaller formal uncertainty than the corresponding line 
pairs, but generally the uncertainties are of the same size 
since P(5/H/3 has a somewhat smaller separation. 

What changes to the derived physical properties of 
the H II regions does the choice of extinction model make? 
This will depend on the wavelengths and strengths of the 
relevant lines, but since the effective optical depths do not 
differ very much among the models it is not expected to 
see large effects on derived physical parameters. This can 
be illustrated by the oxygen abundance calibrator i ?23 = 
([On] -I- [Oiil])/H/3 widely used to determine the metal- 
licity of spiral galaxies. For this ratio the main difference 
will be on the [O ii] lines due to the scattering included in 
the WTC model, while [O ill] only will experience minute 
changes relative to H/3. 

From the uncorrected line fluxes of NGC 604 by Diaz 
et al. (1987) we have calculated the oxygen abundance us¬ 
ing the calibration formula given in Zaritsky et al. (1994) 
for the three cases of extinction models. Extinction cor¬ 
rection assuming a foreground screen and a homogeneous 
slab both give 12-|-log(0/H) = 8.80, while the WTC model 
results in 8.82, i.e. virtually no change for the small dust 
extinction of NGC 604. The spectrum of NGC 5461 is af¬ 
fected by a larger dust extinction and from the data by 
Rayo et al. (1982) we find a difference in oxygen abun¬ 
dance of 0.08 dex between the foreground screen and the 
’dusty nucleus’ model. This is almost negligible compared 
to the estimated precision of the oxygen abundance cali¬ 
bration, roughly 0.2 dex (Pagel et al. 1980), although ex¬ 
tinction correction by the WTC model will systematically 


infer higher oxygen abundances. This result makes the 
widespread use of a simple foreground screen for deredden¬ 
ing justifiable at this magnitude of extinction at least for 
lines not to far apart. For diagnostic lines with longer sep¬ 
arations larger changes are expected, for instance Bautista 
et al. (1995) find deviations of up to 12 % in nebular abun¬ 
dances in a revaluation of Orion Nebula data with an em¬ 
pirical effective nebular extinction law. 

4.J. Spatial variation of the extinction 

Two of the Hii regions, NGC 5461 and NGC 604, has 
sufficiently high signal and large extent that the spatial 
variation of the extinction across an H ii region can be 
studied in greater detail. The E-W aperture (Jul 31) of 
NGC 5461 has been divided into 5 subapertures of length 
3'.'2 corresponding to 84 pc at the distance of NGC 5457 
and the extracted spectrum from each of them treated in 
the usual way. In Fig. 5 the optical depths derived in the 
large aperture and each of the subapertures assuming the 
foreground geometry are depicted at their positions along 
the slit on top of the observed flux in one of the Paschen 
lines, which is less affected by the extinction. As expected, 
the value of ry for the total aperture is close to those found 
in the central subapertures containing most of the flux 
with uncertainties increasing outwards. A scatter in ry of 
~ 0.3 across the region is evident, which could indicate 
a patchiness of the obscuring dust integrated along the 
line of sight of this order. A similar behaviour is found by 
Greve et al. (1989) for the Galactic Orion Nebula. 

We also note that ry increases towards the edge of the 
region (subaperture and 4). This can be interpreted as 
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Spatial extent along slit (aresec) 


Fig. 5. Spatially resolved optical depth across NGC 5461. The 
solid line represents observed flux in the spatial direction along 
the slit in ADU’s on the left ordinate at the central wavelength 
of the P9 emission line. The * denotes the optical depth given 
on the right ordinate with associated vertical error bar as de¬ 
rived from the total integrated spectrum. The size of the ex¬ 
traction aperture is shown as the horizontal dotted line. The 
o’s indicate the optical depths found in the 5 subapertures ex¬ 
tracted over the short dotted lines 


a larger amount of dust residing just outside the brightest 
part of the H ii region possibly due to destruction and out¬ 
flow of dust in the inner starburst environment. It also fits 
well with the picture upon which the WTC dusty nucleus 
simulations are based. Given the uncertainty in modelling 
the dust geometry discussed above one should, however, 
be cautious of drawing any firm conclusion on the detailed 
distribution of the dust. 

The separation of the three bright knots observed in 
the NGC 604 complex are roughly 15"= 50 pc. From the 
overlying screen model we find here somewhat smaller 
variations of ry up to 0.1, while Diaz et al. (1987) also 
found variations in ry of the order 0.3 among 5 areas 
studied in detail. Thus there are signs of variation in the 
optical depths towards giant extragalactic H ii regions on 
scales smaller than 50 pc most likely attributed to a patchy 
distribution of the obscuring dust. 


population synthesis models to the energy distribution 
while taking the dust attenuation into account. 

A substantial scatter between the different sources is 
seen, amounting to more than one unit in optical depths 
for NGC 5461. Part of this can possibly be explained by 
the different aperture sizes and positions used if the dust 
distribution is not totally homogeneous and entirely in the 
foreground, but it may also be related to the problem of 
deriving the extinction from a pair of lines with only small 
wavelength separation. The newer results seem to favour 
the generally high dust extinctions that we derive, but also 
unexplainable wildly deviating values exist. 

Our data display the same tendency as seen in Paper I 
and also commented by Skillman & Israel (1988), that the 
derived extinctions increase with increasing wavelength 
separation. A similar behaviour is found by Bautista et al. 
(1995) in an analysis of the Orion Nebula, where Ay esti¬ 
mated using the lines of a given recombination line series 
(Balmer to Brackett) increases with the mean wavelength 
of that line series. This finding is even based on a value of 
Ry = 5.5, which is appreciable altered from the value of 
3.1 valid for the interstellar medium. (With this value the 
differences would have been more distinct). 

This trend correlates with the findings by Viallefond 
& Goss (1986 and Skillman & Israel (1988) of yet higher 
values of extinction for NGC 598 and NGC 5457, respec¬ 
tively, when comparing radio continuum emission and H/3 
line emission. Both attribute this ’excess’ extinction to 
part of the dust being non-uniformly distributed and/or 
mixed with the gas. Skillman & Israel (1988) point out 
that most of the difference appears to be explainable by 
dust associated with the Hii regions, but there is no un¬ 
ambiguous method to distinguish between the two distri¬ 
butions. It is interesting to see that the absorption at H/3 
for NGC 595 and NGC 5461 from these two studies when 
converted to optical depth in V (ry = 1.4 and 1.8, respec¬ 
tively) match the effective ry of the WTC ’dusty model’ 
in our investigation. 

Similarly Caplan & Deharveng (1986) report and av¬ 
erage excess in Ay of 0.3 mag for regions in the LMC and 
attribute this to clumped interstellar dust and scattering 
dust in a slab near the periphery of the emitting gas. On 
the other hand Caplan et al. (1996) find in a recent study 
of the SMC only a few regions with excess extinction and 
the reddening of the large majority of Hii regions being 
consistent with uniform interstellar extinction. 


4-2. Comparison with previous investigations 

The brightest Hii regions in NGC 598 and NGC 5457 
have been the objective of many extensive studies over 
the years. In the lower part of Table 3 and 4 independent 
extinction determinations available in the literature are 
compiled. The majority are based on the Balmer decre¬ 
ment (Ha/H/3), but for NGC 5461 two other sources are 
included. One is derived from the ratio of Py/H/?, while 
Rosa & Benvenuti (1994) use a different approach fitting 


5. Conclusions 

The observations of giant H ii regions in the spiral galax¬ 
ies NGC 598 of NGC 5457 have demonstrated how the 
dust extinction can be derived on the basis of multiple 
emission lines from the Balmer and Paschen series in the 
optical and near-IR regimes obtained simultaneously with 
ordinary CCD detectors. Although some of these lines can 
only be observed with larger flux errors the wide wave- 
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length interval spanned by the lines and the utilisation 
of the redundancy of the Paschen lines at AA8500-9000 
A leads to a small overall uncertainty compared to the 
use of the two strongest Balmer lines. 

Inclusion of lines beyond 1 /im offers a chance to make 
simple test of the geometrical distribution of emitting 
source and dust in the light of the common but question¬ 
able assumption of a foreground screen. A slab of homoge¬ 
neous mixture can only match the data points for the H ii 
regions with relative low dust extinction. Generally we can 
not discriminate between the simple foreground extinction 
screen and the WTC ’dusty nucleus’ model, although signs 
of deviation from the pure foreground screen are seen at 
near-IR wavelengths. This is in agreement with Calzetti 
et al. (1996), who found that the reddening towards star- 
burst regions in 13 galaxies by large can be explained by 
foreground dust. The choice of extinction model for red¬ 
dening correction of diagnostic line ratios only has small 
effects on the derived physical parameters when lines are 
not widely separated. 

It is, however, clearly seen how the foreground screen 
always yields a lower limit to the dust content and that 
considerable larger amounts of dust can be present in 
the ’dusty nucleus’ configuration. A more definitive dis¬ 
crimination between various geometries, possibly includ¬ 
ing patchiness, and assessment of the actual dust content 
will need investigations at even longer wavelengths in com¬ 
bination with optical lines (Puxley & Brand 1994; Genzel 
et al. 1995; Calzetti et al.l996) and more refined mod¬ 
elling (Witt & Gordon 1996). While a lot can be learned 
from high-S/N observations of Py obtainable with the 
high quantum efficiency of the new generation of CCD’s, 
it is desirable to have data on emission lines further out in 
the near-IR, e.g. P/3, Bry, Bra, which can not be observed 
with the same instruments as the blue optical lines. The 
kind of studies presented in this paper can serve to bridge 
the two wavelength ranges. 

One obvious result from the present study is the 
demonstration of a spatial variation in the optical depth 
towards the emission regions of the order 0.2-0.3 on small 
angular scales. More detailed investigations of one of the 
HII regions suggest that the extinction has it highest value 
at or close to the edge of the region. This can explain 
at least some of the disagreement in the quoted extinc¬ 
tion of HII regions by different authors. It also emphasizes 
the importance of sampling exactly the same slit position 
and aperture size if combining data from different obser¬ 
vations, and the care that most be taken when comparing 
physical parameters from miscellaneous sources in the lit¬ 
erature. 
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